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Abstract. The level spacing distribution is numerically calculated at the disorder-induced 
metal-insulator transition for dimensionality d = 4 by applying the Lanczos diagonalisation. 
The critical level statistics are shown to deviate stronger from the result of the random matrix 
theory compared to those of d — 3 and to become closer to the Poisson limit of uncorrelated 
spectra. Using the finite size scaling analysis for the probability distribution Q„(E) of 
having n levels in a given energy interval E we find the critical disorder W c = 34.5 ±0.5, the 
correlation length exponent v = 1.1 ± 0.2 and the critical spectral compressibility k c » 0.5. 
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1 Introduction 

The statistical description of electronic spectra at the Anderson transition is one of 
the subjects of central interest. This quantum phase transition implies that there 
exists a non-vanishing value of the disorder of the random potential, at which the 
system undergoes the crossover between a conducting and an insulating phases with 
delocalised and localised electron states, respectively. The spatial dimensionality d 
of the system plays an important role in determining the critical properties of both 
the conductivity and the level statistics. According to the one-parameter scaling 
theory of localisation 0] all single-electron states in the one and two dimensions 
are localised even for arbitrarily weak randomness, provided that the time-reversal 
and the spin-rotational symmetries are preserved. Therefore, the lowest integer di- 
mension for which the disorder-induced metal-insulator transition (MIT) occurs for 
non- interacting particles is d = 3. Numerous computer simulations performed on 
transport electron properties have confirmed the scaling hypothesis (see, for exam- 
ple, [|| U and references therein). Furthermore, the MIT has been also found at the 
next higher integer dimensionality d = 4, by applying the transfer-matrix method 
(TMM)§,§. 

The energy level statistics in two dimensions (2D) do not exhibit critical behaviour || 
f?J , while in 3D they do || || |l0| in agreement with the earlier results obtained by the 
TMM H . A new universal scale-independent level spacing distribution P(s) has been 
predicted exactly at the MIT in 3D Q . This distribution differs characteristically from 
the results of the 'classical' random matrix theory (RMT) by Wigner and Dyson [IT], 
|i~2"[ , which is valid for weakly disordered conductors [jl3], Q . An analytical approach 
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has been developed recently |l5|, [l6|, which supposes that the shape of the P(s) at 
the mobility edge is given by a combination of the dimensionality d and the critical 
exponent v of the correlation length. At present this relation is intensively investigated 
for models of different basic symmetry, namely for the orthogonal jl?], [18], |l9), the 
unitary |2(], and the symplectic (2^, ^3], [24] critical ensembles. Here the critical 
disorders and the exponents are known, and the dependence of Pis) on v and d could 
be determined quantitatively. In order to obtain further insight, it is of great interest 
to study also critical parameters and the level statistics at higher dimensions. 
In this paper we report results of detailed numerical calculations on the level statistics 
using the Anderson model for the 4D hypercubic lattice. We show that the nearest- 
neighbour level spacing distribution P(s) exhibits critical behaviour. Its shape is size- 
invariant at the transition, and its form considerably differs from that in 3D. Using 
the level statistics method in terms of the complete distribution of having n levels in 
a given energy interval we detect the metal-insulator transition, which corresponds 
to the disorder W c = 34.5 ± 0.5. We perform the finite size scaling analysis in order 
to find the one-parameter scaling function and to determine the correlation length 
exponent v, which is smaller than that in 3D. Combining the results obtained for 
the statistic Jo (defined below) and the variance of the level number we argue that 
the critical spectral fluctuations are stronger in 4D than in 3D, being closer to the 
Poissonian statistics. 

2 Model and computational procedure 

The model for the 4D disordered system is defined by the Anderson Hamiltonian 

H = Y1 e n C l C n + V ^ ( C n c m + C„4j, (1) 
n n^m 

where (c„) is the creation (annihilation) operator of an electron at a lattice site n, 
and m denotes the sites adjacent to the site n (their number equals 8 for the simple 
hypercubic lattice). The site energies e„ are randomly distributed according to a box 
distribution with a width W, which plays the role of the disorder parameter. The 
second term describes the hopping between the nearest-neighbour sites in the lattice. 
Our considerations are restricted to the particles without spin and with no magnetic 
field. It was earlier found by the TMM Q that the MIT at the band centre E=0 is 
close to the disorder W = 33.2. For 3D the critical disorder is smaller, W « 16.5 [Q. 
After numerical diagonalisation of the Hamiltonian (|l|) with periodic boundary con- 
ditions using the Lanczos algorithm we obtained the exact discrete spectrum of the 
electrons for simple 4D hypercubic lattices of various sizes ranging from L 4 = 4 4 to 10 4 
and for different disorders W . Linear and energy scales are measured in units of the 
lattice constant and the overlap integral between adjacent sites (V = 1), respectively. 
The levels were taken from energy intervals centred at E — so that they belong to 
the critical energy region, defined by the condition L < £ = (\E — E c \/E c )~ u , where £ 
is the correlation length and E c is the mobility edge. The number of realizations for 
a given size L was such that the total number of eigenvalues amounted as much as 
10 5 . We have checked that the density of states p — (AL 4 )^ 1 around the band centre 
slightly varies with the energy (A is the mean level spacing). Therefore, the careful 
unfolding procedure has been applied for the spectra of all pairs of {W, L}. 



2 




Fig. 1: The nearest neighbour level spacing distribution P c (s) for d = 4 at the 
disorder W = 34.5 for different lattice sizes L 4 . Full line is P c (s) for d = 3 taken 
from |19|. Dashed and dotted curves are Pgoe(s) (§) and Pp(s) (|^), the distri- 
butions for the metallic and insulating limits, respectively. 



3 The critical level spacing distribution P c (s) 

A traditional way to study the statistical properties of disordered spectra is to consider 
the level spacing distribution P(s) which is defined as the probability density of 
nearest neighbouring levels. It is known that in the metallic region P(s) is very close 
to the Wigner surmise |b| for the Gaussian orthogonal ensemble (GOE) of random 
matrices [ pd| , namely 

^goe(s) = |.scx P (-J S 2 ), (2) 

where s is measured in units of the mean level spacing A. In the localised region 
the energy levels are completely uncorrelated, and hence the spacings are distributed 
according to the Poisson law 

Pp{s) =exp(-s). (3) 
In a similar way as for the 3D Anderson model || ||, 10 1, the third, universal distri- 



bution P c (s) is supposed to be revealed exactly at the critical point, which is different 
from both of the above laws. For systems of finite size the level statistics is expected 
to change continuously from Pgoe(s) through P c (s) to Pp(s), when increasing the 
strength of fluctuations of the random impurity potential. When increasing the size of 
the system the distribution tends towards either Pgoe(s) (@) or Pp(s) (||), depending 
on whether the disorder is below or above its critical value, respectively. 
Fig. H shows the function P(s) calculated at the disorder W very close to the critical 
point. The fact that the data within numerical errorbars lie on a common curve inde- 
pendent on L, defines the critical level spacing distribution P c (s) (see also section 5). 
In comparison with d = 3 the critical P c (s) for d = 4 is closer to Pp(s). It is worth 
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Fig. 2: Large-s part of the critical level spacing distribution P c (s) of the Fig. 1 for 
various system sizes L. Dashed-dotted straight line, In P c (s) = —1.4s, is the best 
fit to the data. Dashed and dotted curves are Pgoe(s) and Pp(s), respectively. 



noticing that our numerical data deviate considerably from the interpolation formula 
P c (s) = B s exp(— A s 1 + 1 / du ) derived analytically in |l6| ], where the coefficients A and 
B are defined by the normalisation conditions. 

For small spacings we found the linear increase P c (s) = Bs, as consistent with the 
orthogonal symmetry. However, the linear slope B « 2.1_Bgoe is larger than in 
d — 3 H and, consequently, in the RMT. This indicates that the level repulsion 
becomes weaker with increasing d. The behaviour of P c (s) at larger spacings is well 
described by the sub-Poissonian form P c (s) oc exp(— As) with A = 1.4 ± 0.1. One 
can see in Fig. || that P c (s) changes by several orders of magnitude in the interval 
2 < s < 10. As expected, this asymptotic decay is slower compared to d — 3, where 
the exponential rate is A « 1.9 (l9), but faster than the Poissonian decay. 



4 The probability distribution Q n (s) 

In order to perform the finite size scaling analysis and to study how the level statistics 
behave around the critical point, we have calculated the dependence of the statistics 
of neighbouring spacing on the disorder W for different L. To include the entire range 
of spacings, we use here the probability distribution of having n eigenvalues in a given 
energy interval of the width s: 

/>oo />oo />oo 

Q n (s)= I n (s')ds'= ds' p n (s"),ds", (4) 



where I n (s) being the cumulative distribution of n successive levels. The function 
Q n (s) is known for the metallic limit from the RMT |l2|]. In the insulating limit the 
Poisson process Q n (s) = s"exp(— s)/n\ governs the completely uncorrelated spec- 
trum. 

In what follows we investigate the probability to have no level, n = 0, within the 
bin s, which implies the level statistics for the case of nearest neighbour spacings. In 
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Fig. 3: The probability distribution AQo(s) at different disorders W = 20 (upper 
set), 33.2 (middle set) and 40 (lower set) for various sizes L of 4D hyper-cubic 
lattice. Solid line is the RMT result taken from the Table A-19 of Ref. Jl2| . The 
arrows indicate the direction of the size effect on either side of the transition. 



order to show (L, VF)-dependence we plot in Fig. g the absolute deviation from the 
Poisson process AQ ( S ) = — [Qo{s) — exp(— s)]. The data for the metallic phase, i.e 
for W < W c , turn out to be closer to the GOE. The larger the system size, the closer 
the level statistics to this limit. When W > W c the data approach zero, however with 
the opposite size effect. Closer to the critical disorder all of the data within statistical 
uncertainties start to fall onto a common intermediate curve independent of L. This 
is the manifestation of critical behaviour. Similar scaling properties are also observed 
for Q n>0 (s). 



4.1 The spectral statistic J„ 

To paramctcrisc the distribution Q n (s) we deal with a global statistical quantity 
defined as follows 

poo 

J n = / Q n (s)ds, (5) 
Jo 

which takes the whole range of the spacings into account. For n = one can easily 
show that it is related to the spacing variance (S 2 s) as 

Jo = l(s 2 ) = l((6 2 s) + l). (6) 

For investigating the critical properties of the spectra at the Anderson transition the 
statistic Jo was for the first time introduced in []l7t . It proved to be more efficient than 
those which have previously been used in studying the 2D and 3D cases J(| ||, |To| . 
This is due to that since the probability density P(s) = po(s) and the cumulative 
distribution Iq(s) are normalised to unity, one had to choose some spacing s* in order 
to weigh the functions after or before this point M, N, MO] for a given pair of the 
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Fig. 4: Disorder dependence of Jo for various sizes L. Full line is the RMT 
result J GOE = 0.643 [Q. Inset: the derivative dJ /dW near W c vs. size L. Fit of 
numerical data shown by dashed line yields v = 1.1. For comparison, the slope 
of dotted line corresponds schematically to v — 1.45 for d = 3 [J1C|]. 



parameters {W, L}. This leads however to loosing part of the information obtained 
from the diagonalisation and, as a result, decreases the accuracy. In contrast, the set of 
parameters J„ does not require s*. Therefore, J„ is more successful, in our opinion, 
for demonstrating the scaling properties of the level statistics with less number of 
realizations and, as consequence, to locate the MIT more precisely. For the GOE 
Jo = 0.643, Ji = 0.922, J n=00 = 1 |jl2). In the localised regime, one has simply 
J n = 1 for any n. For d = 3 the set of the critical numbers J^ has been obtained 
in |l7j], for example, Jq w 0.714 and (s 2 ) rj 1.42. In the presence of strong spin- 
orbit interactions, where the MIT occurs even for d = 2, another set describes the 
critical symplectic ensemble ||^] . 



4.2 The finite size scaling and the critical exponent 

We have calculated the dependence of Jo on W near the transition for various system 
sizes. All computed data lie within the interval between Jq oe and Jq , gradually 
growing from the former to the latter limit, when increasing the disorder W. The 
increase of Jo develops faster with L. For an infinite system this change would trans- 
form to a discontinuous crossover between these two limits exactly at the transition 
W c . One observes from Fig. |4|, that Jq{W) exhibits critical behaviour. The common 
crossing point Jq w 0.79 ((s 2 ) ~ 1-57), where the size effect on the statistics changes 
sign, corresponds to the transition. Intersection of the curves of different L enables 
us to determine the fixed point more precisely, W = W c = 34.5 ± 0.5. The obtained 
value is in reasonable agreement with that computed previously by the TMM Q H . 
However it markedly deviates from the linear relation W c (d) = (d — 2)W c (d = 3) Qj. 
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Fig. 5: One-parameter scaling dependence of the statistic Jo on L/£ for different 
system sizes L and disorders W. Full straight line is the RMT result Jo = 
0.643 O. Dotted line is (Sh with J C = 0.79 and the coefficient C = 0.043. 



At the fixed point the correlation length diverges with the exponent v. 

£(W) oc \W- W c \- V . (7) 
By using the lowest terms of the expansion for the function Jo (W, L) near W c 

J (W, L) « J C + C(W - W C )L^ V = J C +C (jj ' , (8) 

one extracts the critical exponent of the correlation length. Inset of Fig. || shows that 
the data within the numerical errors are well described by the linear approximation 

i dJ (W,L) _ x 

ln dw lnL - (9) 

The estimated value v = 1.1 ± 0.2 is consistent with previous findings It appeared 
to be smaller than that for d = 3 [y w 1.45 0|). On the other hand, it is still larger 
than the standard mean- field result %p = 1/2, valid for the upper bound d u of the 
Anderson transition ^BJ. Based on our results one can argue that d — 4 is definitely 
lower than d u , which is believed to be equal infinity [^6| . 

By introducing the scaling variable £, which is identified as the correlation length, 
it is possible to replot all of the data in Fig. || into a single-parameter function, as 
shown in Fig. The resulting curve Jo(L/£) consists of two branches characteristic 
of different regimes of the disorder-induced MIT. The decaying branch corresponds 
to the metallic phase (W < W c ) and the growing one belongs to the insulating phase 
(W > W c )- Outside of the critical region, i.e. L/£ > 1, the numerical data deviate 
from the linear approximation (|J). The one-parameter finite size scaling procedure 
allows one to find the dependence of the correlation length £ on the disorder. Close 
to W c the numerical results of £(W) shown in Fig. || give a satisfactory agreement 
with (0), while far apart from the fixed point one observes a considerable discrepancy. 
Thus, the critical behaviour of the spectral statistics in 4D is typical for the Anderson 
transition and analogous to that for the 3D case || |[ |l0) . 
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Fig. 6: Disorder dependence of the correlation length £(W) for 4D disordered 
system. Dotted line is (Q). 



5 Variance of the energy level number 

While the spacing distribution P(s) and, consequently, Jo probe the short-range cor- 
relations in the electron spectra, the variance of the number of levels (S 2 N(E)) in a 
given energy interval E, which describes the global spectral rigidity, can provide the 
information about fluctuations on scales much larger than A. It is defined as a width 
of the distribution of N levels in the interval E: (S 2 N(E)) = J2n=o( n ~ ( n ))Qn(E). 
In the extreme insulating limit W 3> W c the number variance obeys the ordinary 
Poisson law (S 2 N(E)} — (N(E)). For any finite W the eigenfunctions can spatially 
overlap, so that the fluctuations (S 2 N) are reduced below the Poisson limit due to 
level repulsion. When the system is a good conductor, i.e W <C W c , the electron states 
are spreaded over the entire volume. Therefore the number variance for (N) 3> 1 can 
be approximated by the Dyson formula 

(5 2 N) = A in(jv) + 7 , 7 « 0.44, (10) 

valid for the Gaussian orthogonal ensemble of random matrices ]l2, |^|. Thus, the 
relative variance (S 2 N)/(N) for large (N), which is also known as a spectral com- 
pressibility k, changes from zero to unity as W increases. Of particular interest is the 
question how does the relative variance behave at W = W c . For instance, in 3D sys- 
tems it was numerically shown that its critical value equals a constant k c rj 0.27 p7| . 
It may also be worth considering in more detail the number variance for weakly dis- 
ordered 4D systems in order to compare with the results of the diffusive theory [JL4J . 
Fig. [?] shows the numerical results of the relative number variance around the transi- 
tion. The size dependence below and above the critical point is very similar to that 
of AQo(s) in Fig. ||. With increasing L the data approach the Dyson result jic| ) 
if W < W c , while k — > 1 if W > W Cl respectively. At W = W c = 34.5 the ratio 
(S 2 N)/(N) is almost insensitive to L. This justifies again the existence of the in- 
termediate scale-invariant statistics at the long-range energy correlations. The ratio 
(S 2 N)/(N) at the MIT decreases very slowly with increasing the mean level num- 
ber (N). 



8 



1 




o 1 1 — 

1 10 

<N> 



Fig. 7: The relative number variance (5 2 N)/(N) of the spectrum as a function 
of the mean number of levels in a given energy interval (N(E)} for various sizes 
L at three disorder values: W — 20 (lower set), 34.5 (central set) and 47 (upper 
set). Full line corresponds to (tlfl) . 



Recently the analytical theory has been proposed |27| suggesting that the compress- 
ibility at the mobility edge is totally characterised by properties of the critical eigen- 
states and is directly determined in terms of the multifractal exponent /i = d — D2 

d{S 2 N) 11 1 
n c = hm —, — = — < -. (11) 

The results of \x and k c obtained numerically so far for lower dimensions d < 3 satisfy 
reasonably this relation. Our 4D-data for sufficiently large values (N) > 20 yet obey 
at least the condition for the above upper limit, going down to k c w 0.45 — 0.5. This 
value appears to be larger than that of 3D, indicating that the critical spectral rigidity 
diminishes. However the accuracy of (S 2 N) is still not high enough to reach a precise 
saturation value at (N) 1 and, as a result, to provide a reliable estimate of fi 
from (|ll]). In addition, an independent computational analysis of the multifractality 
of eigenfunctions in 4D would be needed. 



6 Conclusions 

We have numerically calculated the critical distribution of neighbouring spacings P(s) 
at the metal-insulator transition for the 4D Anderson model. The disorder-induced 
crossover between the Wigner and the Poisson statistics for finite systems is shown 
to obey a one-parameter scaling law as for lower dimensions. The finite size scaling 
analysis of the level statistics allows one to locate the critical disorder W c and to 
detect the correlation length exponent v. Comparing the obtained results with lower 
dimensions we are led to the conclusion that the spectral correlations at criticality 
depends on the spatial dimensionality, becoming weaker with increasing d. The fur- 
ther systematic study of the dimensionality dependence would be desirable in order to 
answer the question how the critical level statistics favours the Poissonian limit, when 
approaching the upper bound of the Anderson transition |2(|. Another interesting 
problem for simulations in 4D could be spectral fluctuations in a weakly disordered 
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metal to check the non-perturbative theory 
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